Ideal glass transitions for hard ellipsoids. 
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For hard ellipsoids of revolution we calculate the phase diagram for the idealized glass transition. 
Our equations cover the glass physics in the full phase space, for all packing fractions and all aspect 
ratios Xq. With increasing aspect ratio we find the idealized glass transition to become primarily be 
driven by orientational degrees of freedom. For needle or plate like systems the transition is strongly 
infiuenced by a precursor of a nematic instability. We obtain three types of glass transition lines. 
The first one {(f>c ) corresponds to the conventional glass transition for spherical particles which is 
driven by the cage effect. At the second one {4>c ) which occurs for rather non-spherical particles 
a glass phase is formed which consists of domains. Within each domain there is a nematic order 
where the center of mass motion is quasi-ergodic, whereas the inter-domain orientations build an 
orientational glass. The third glass transition line {<l>c ) occurs for nearly spherical ellipsoids where 
the orientational degrees of freedom with odd parity, e.g. 180° flips, freeze independently from the 
positions. 
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I. INTRODUCTION 

The dynamics of a molecular system which is super- 
cooled towards the glass transition shows a variety of phe- 
nomena related to the nontrivial interplay between ori- 
entational and translational degrees of freedom caused, 
e.g. by steric hindrance. 

Already in thermodynamic equilibrium molecular sys- 
tems show, compared to simple liquids, a variety of differ- 
ent physical behavior. At low enough densities (or high 
enough temperatures) they form an isotropic liquid. On 
increasing the density they can undergo a transition into 
a crystal or several different liquid crystalline phases (like 
e.g. a nematic phase). A crucial part of the interaction 
which causes these phenomena is given by the shape of 
the molecules. Therefore one may also expect new char- 
acteristic features for the glass transition in such systems. 
A model system which allows to study the translation- 
orientation interplay is a system of N hard ellipsoids of 
revolution in a box of volume V. The fluid of ellipsoids 
is characterized by two parameters: The aspect ratio 
Xq = - relating b and a, the major and minor axis of the 
ellipsoids, and the packing fraction (j> which is related to 
the number density p = y hy (p = ^nXop. 

In this work we start from a theory of liquids. We 
demonstrate that a single set of equations allows to de- 
scribe the glassy behavior for almost spherical particle up 
to large aspect ratios where pre-nematic order becomes 
crucial. 

The choice of the model system has also been moti- 
vated by the successful application of the ideal mode cou- 
pling theory (MCT) for simple liquids ^ , particularly to 
neutral colloidal suspensions. MCT gives a closed set of 
equations for the intermediate scattering function S'(q, t). 
Comparison between experimental H and MCT-results 
|Q,|| have shown good agreement for the colloidal sys- 
tems which usually are modeled by hard spheres. Fur- 
ther tests of the MCT for other systems can be found in 



e.g. refs. y-|§l. MCT in its original form describes an 
idealized glass transition which is indicated by breaking 
of ergodicity, at a critical density p^ (or critical temper- 
ature Tc). The corresponding non-ergodicity parameter 
/(q) — hmj^oo 'S'(q, i)/S'(q, 0) becomes nonzero at Pc (or 
Tc). Recently the mode coupling equations have been 
extended to molecular systems. The dynamics of liquids 
of rigid molecules composed of M atoms can either be 
described by site-site correlators 5q,,/3,q;,/9 = 1,2, ..., M 
or molecular correlation functions Simn,i'm'n' (q, t) where 
for the latter one decomposes the degrees of freedom into 
the center of mass and orientational components (see e.g. 
[0,||). The density p{x,fl,t) is a function of the center 
of mass coordinate x and the orientation il — ($,6',x) 
which is specified by the three Eulcr angles. Expanding 
/9(x, Q,, t) with respect to a product basis of plane waves, 
giqx ^j^^ generalized spherical harmonics D5„„(J7) one ar- 
rives at the tensorial density pimn (q, t) . I runs over all 
positive integers including zero and m and n takes integer 
values between —I and I. Then the molecular correlators 
are defined as follows : 

Slmn,l'm'n'{(l,t) = — (p;,„„ (q, t)p;/,„/„/ (q, 0)) (1) 

The extension of MCT to molecular systems has been 
done for the molecular representation for a single dumb- 
bell in a simple, isotropic liquid ^, for molecular liq- 
uids of linear molecules pO| and for arbitrarily shaped 
molecules by use of nonlinear fluctuating hydrodynamics 
pT| and by Mori-Zwanzig projection formalism |T2]. An 
MCT-approach using a site-site description has recently 
been worked out |l3]. Because a hard ellipsoid corre- 
sponds to a rigid body with infinitely many constituents 
it is the molecular representation which is the only ap- 
propriate one. Since we consider ellipsoids of revolution 
the third Euler angle x becomes redundant. This means 
that we have to consider Sim,o,i'm'o{^'i t)i only. Using the 
q-frame |8|, i.e. one chooses q = (0,0,(7) = Qo where 



q— |q|, these correlators become real and diagonal in m 
and m' |lo|: 



S, 



ImO.l'jn 



o(q,0 



'Sw{q,m,t) 



(2) 



The head-tail symmetry of the ellipsoids implies that 
these correlators vanish for I + I' odd. For given Xq the 
critical packing fraction (pcO^o) can be determined by cal- 
culating the (unnormalized) non-ergodicity parameters 
Fii> {q, m) = limt_oo 5";;' (g, m; t). 



II. MOLECULAR MODE-COUPLING 
EQUATIONS 

Using the densities pim{q,t) and longitudinal transla- 
tional currents jf^{q,t) and rotational currents Ji^{q,t) 
as slow variables set for the Mori-Zwanzig projection 
operator technique the molecular MCT-equations have 
been derived and can be found in refs. [Eo| and ||l^. The 
time dependent molecular MCT-equations can be repre- 
sented as follows: 

-S((7, m, t) = N^(g, m, t) + N^{q, m, t) 

^N"{q,m,t) = ~nl{q,m)Siq,m,t) 

- ^ i^aa' {q, m)N" {q, m, t) 

a.' 

-02(g,m)V f m"'"'(q,m,i-i')N"'(q,m,i)di' 



and: 



m^]? {q, m, t) 



qi-q2 mi ^1,^2 



/-m.irix ,a,Q 



{q,qi,<i2) 



1^ 1^ y^ 

mi !i,!2 ''''i-'b 

'i''2 

Sh,i[{qi,mi,t) Si.^M^{q2,m2,t) (4) 



The indices a, a' G {T, R} refer to either translational 
or orientational currents while Ot and Or are the mass 
and moment of inertia, respectively. N" (g, m, t) are the 
current-density correlation functions for translational 
(a = T) and rotational {a = R) currents multiplied 
with q and ^yl{^+l), respectively. The microscopic fre- 
quency matrix is denoted by fta{q, ra) and is determined 
by the static molecular correlators. In the absence of 
memory effects (m"'" = 0) the equations are just a set 
of coupled harmonic oscillators with friction Vaa' for vi- 
brational (a — T) and rotational {a = R) oscillations. 
E.g. the translational mode with I = 1' = is the propa- 
gating phonon mode and the modes with I = I' > which 
exhibit a frequency gap at g = are localized oscillators. 
For m"'" ^ nonlinearities occur. Their physical ori- 
gin are memory effects. The corresponding memory ker- 
nel is a correlation function of fluctuating forces. Since 
fluctuating forces can decay into a pair of density excita- 
tions this kernel is approximated as a sum of all possible 



bilinear products of density correlation functions. Such a 
non-linear feed back mechanism can cause an ideal glass 
transition with non-trivial dynamics. The glass transi- 
tion for eqs. (J^) and (|j) is investigated in the following 
part of the paper. The explicit expressions for the ver- 
tices V"" for arbitrary q can be found in ref. ||l^ and for 
the q- frame in ref. jlj]. The vertices V"" only depend 
on the static correlators Sw (g, m) and the direct corre- 
lation function cii> (g, m) which are related to each other 
by the Ornstein-Zernike equation. We have determined 
cii'{q,m) within Percus-Yevick approximation. 

It has been shown that the liquid phase of hard ellip- 
soids is well described by these approximations |14|. Al- 
though Percus-Yevick theory fails to describe crystalliza- 
tion it yields a nematic instability |l4| which is in reason- 
able agreement with Monte-Carlo simulations |l^ even if 
Percus-Yevick theory still underestimates the tendency 
towards orientational order. This nematic instability will 
play an important role in the following. For the solution 
of the Percus-Yevick equations we have chosen a cut-off 
Ico — 4 for I and I', whereas the MCT equations were 
truncated at Ico = 2. We are confident that even such a 
small number of molecular correlators enables us to cap- 
ture the correct physics of the transition. 



III. SOLUTION OF THE MMCT FOR HARD 
ELLIPSOIDS 



(3) 



The numerical solution of the eqs. (p|)-(H) for t — > 
oo yields the non-ergodicity parameters Fu' (g, m) — 

F(g, m). In the limit of t > cx) the following set of 

nonlinear equations for F(g, m) has to be solved in an 
iterative way: 

^gr(g)g?'(g) (.F(g,m)-y" S-\q,m)F{q,m) 

aa.' 

+ F(g,m)-S(g,m) = (5) 

qi.q2 mi ii,i2 '''''r'2 
' ''r'2 

Fh,i[{qi,'mi) Fi^^i'^{q2,m2) (6) 



with 



irii) 



,/iii + i) 



for a = T 
for a = R 



(7) 



with ^{q, Tn) = 

lirriz-^o — zm"'" (g, to, z) = limt^^m"''" (g, to,, t) we de- 
note the long time limit of the memory kernel. From a so- 
lution of these equations we obtain the phase diagram for 
ideal glass transitions which is shown in fig. |l|. This figure 
also contains two dashed-dotted lines ^nemi^o) indicat- 
ing the location of the nematic instability as it arises in 



thermodynamic equilibrium from PY-theory |L4[ . These 
two lines are in agreement with density functional the- 
ory |1^ and Monte Carlo simulations [|5|. In addition 
there are three glass transition lines, each for Xq < 1 
and Xo > 1. First of all we will discuss the critical line 
(f>c (-^^o) (thick solid line) at which both translational 
and orientational degrees of freedom for / and /' even 
undergo a discontinuous ergodic to non-ergodic transi- 
tion (also called type-B transition). The existence of 

(t>c (Xq) has been established for 0.35 < Xq < 2.5. In 



detected. 

Probably the most interesting result is the third criti- 

( Ft' ^ 

cal line (/>c (-'^o) (dashed line) which is shown schemat- 
ically in figure ^ for Xq > 2.0 (prolate ellipsoids) and 
Xq < 0.5 (oblate ellipsoids). In this region the glass 
transition lines are close to the nematic instability line. 
The existence of (j)c (Xq) is based on our following ob- 
servations. Increasing (p for 2.1 < Xq < 2.5 we find a 

glass transition at 0c (-^o) where all Fui{q,m) become 
nonzero. Increasing (j) further we find in addition a see- 



the region where the (j)[ {Xq) glass transition occurs the ond solution i^/;, {q, m) for > (/)i^' (Xq). This solution 



equilibrium system shows crystallization. Being a first 
order phase transition the onset of crystallization give s 
two densities (e.g. from Monte-Carlo simulations |15|) 
resulting from a Maxwell construction. The 0c (Xq) 
glass transition line is well bracketed between these two 
densities. This indicates that the mode coupling equa- 
tions describe a glass transition in the meta-stable region 
of a super-cooled liquid. The physical origin of the glass 
transition depends strongly on the location on (j^c (-'^o)- 
For aspect ratios Xq close to one the transition is domi- 
nated by the center of mass correlator SqqIq, 0). 

To illustrate this point we have plotted in fig. 
the static center of mass correlator and the "quadrupo- 
lar" correlator 6*22(9,0) and their corresponding non- 
ergodicity parameters for Xq = 1.3. This was done di- 
rectly above the critical packing fraction 4>c — 0.549. The 
first peak at qmax = 9 ~ 6.6[a~^] of S'oo(<Z, 0) dominates 
the transition which is the manifestation of the cage ef- 
fect. Stronger deviations of the ellipsoids from spherical 
symmetry, however, alter this behaviour. This is demon- 
strated in fig. ^ where we have plotted the same correla- 
tors as in fig. gbut for an aspect ratio of Xq = 2.3. Now 
the peak at g « of the quadrupolar correlator S22{q, 0) 
(which is for q=0 the Kerr constant for non polar fluids) 
dominates the breaking of ergodicity. The half width 
Ag (at half maximum) of this peak defines a correlation 
length ^ = 2n/Aq. In fig. we have plotted ^ at the 

glass transition line (either 0c or 0c ) as a function 
of the aspect ratio Xq for prolate ellipsoids. 

Within the glassy phase, i.e. for (f) > (j)}, [Xq) a con- 
tinuous (also called type-A) glass transition occurs at 
the critical lines 0c (-^0) (thin solid lines) at which the 
self part of the correlators with I and V odd freezes. This 
type-A transition can only occur if the corresponding 
vertex V is large enough. For this to happen the aspect 
ratio should clearly be different from one. The reader 
should note that four points (open circles) were deter- 
mined exactly and the thin solid line is schematic, show- 
ing that 0c {Xq) has to increase if Xq is changed towards 
one, in order to keep the vertices large enough. The phys- 
ical interpretation is that at 0c {Xq) the 180°-jumps of 
the ellipsoids become frozen. This resembles the forma- 
tion of orientational glasses. One possible candidate for 
such a transition might be plastic crystals like the carbo- 
ranes [0| although presently only type B transitions are 



has the feature that i^/^, (g, m) is practically zero with 
the exception of a well pronounced peak for F22 {q, m) at 
q = and with a width of order ^. We have shown this 

in fig. pla where we have plotted F22{q,0) at the 0c 
transition for Xq = 2.5. This is plotted together with 
the static correlator 522(9,0) and the normalized func- 
tion f^2i'1.0) = F^2(g,0)/522(<Z,0). In fig. |b wc have 
plotted the same quantities for Z = T = 0, the center of 
mass correlator. /oo('?, 0) does not exceed 0.15 although 
the orientations are frozen. This means that the system is 
" quasi-ergodic" in the sense that for length scales I <C £, 
the ellipsoids show a (nematic) orientational order and 
the center of mass behaves "quasi-ergodic" decaying to 
a very small value. Whereas for length scales I ^ £, the 
orientations as well as the positions have non-decaying, 
long-time correlations and are frozen. The easiest way 
to think of such a system is due to the formation of liquid 
crystalline (nematic) domains with a size of the order of 
^. This is visualized in fig. g. Fig. ^ shows that the do- 
mains can be quite large. For Xq = 2.5 (the aspect ratio 
where the type-B' transition occurs first and where there- 
fore the type-B' transition with the smallest domain size 
shows up) we obtain from our calculation a domain size 
with ^ « 30. Within the domains the center of mass is 
"quasi-ergodic" i.e. liquid like whereas the orientations 
are frozen with a nematic order. In our idealized MCT 
an ellipsoid can not move from one domain to the other. 
In connection with this it is also interesting to mention 
that two types of type B transitions were also found for 
the center of mass correlator of a simple liquid of hard 
spheres with an attractive interaction given either by the 
Baxter model [|l8[ or a Yukawa potential [|l9[. The ex- 
istence of these two solutions for Xq > 2.1 reflects the 
competition between the frozen positional disorder due 
to the cage effect and the tendency to form a nematic 
phase. Since Fli{q,m) < Fii{q,m) for 2.1 < Xq < 2.5 
the second solution is unphysical [pO|. However, for 

( H') i—r 

Xq > 2.5 and 0c {Xq) < < 0nem(A^o) WC only flnd 
one solution which has all the features of F^, {q, m) de- 
scribed above. We stress that the existence of the criti- 
cal line 0c (^0) trusts on our choice of slow variables 
which includes the nematic order parameter and there- 
fore accounts for the occurrence of the weakly first or- 

der nematic transition. Since 0c {Xq) is rather close 
to 0„em(A'o) "quasi-critical" fluctuations appear which 



also slow down the entropy fluctuations. We do not think 
that those will qualitatively change the phase diagram. 
On the other hand the concept of a glass transition in- 
duced by the vicinity of a second order phase transition 
was already introduced by a MCT approach [pl| in or- 
der to describe the experimentally observed central peak 
phenomenon close to a ferroelectric instability. 



IV. CONCLUSION 

In conclusion, we have shown that hard ellipsoids ex- 
hibit a rather intriguing phase diagram obtained from the 
idealized mode coupling theory for molecular systems, 
where the orientational degrees of freedom and their cou- 
pling to translational ones is incorporated. Particularly, 
we predict a glass transition for Xq > 2 which is driven 
by a precursor of a nematic phase. Ellipsoids show two 
type-B glass transition lines {(pc and 0c '). One, (j)[ , 
is dominated by the cage effect whereas the other one 
(j)c is caused by an orientational (nematic) instabil- 
ity. Besides this a type-A glass transition occurs for al- 
most spherical ellipsoids where the orientational degrees 
of freedom with odd parity, e,g, 180° -flips freeze inde- 
pendently from the translational ones. It would be very 
interesting to check these predictions by experiments or 
simulations. 
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FIG. 1. Phase diagram for the ideal glass transitions. The 



horizontal axis shows Xo scaled with 



x^+i 



The type-B glass 

transition lines 4>c{Xo) and 4>c (Xo) (see text) are depicted as 
thick solid and dashed lines, respectively. The thin solid line 
is the cj>c (Xq) glass transition line. The nematic instabil- 
ity occurs at (jinem{Xo) which is shown as thin dashed-dotted 
lines. The inset shows the situation around Xo — 2.5 where 
the 4>c{Xo) glass transition line merges into the (j)c {Xq) tran- 
sition line. 
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FIG. 2. The static structure factor Sii'{q,m) is plot- 
ted together with the non-ergodicity parameter Fhi {q, m) for 
Xo = 1.3 and = 0.549 (directly above the non-ergodicity 
transition). Part a) shows the center of mass correlator 
I = I' = m = Q whereas part b) shows the quadrupolar corre- 
lator 1 = 1' = 2,m = Q. 



FIG. 3. A similar plot as fig. H is shown but for Xo = 
2.3 and cj> = 0.617 (again directly above the non-ergodicity 
transition) which is already close to the nematic instability. 



FIG. 4. At the glass transition the correlation length for 
parallel orientation obtained from the half-width of the peak 
of S'22(?,0) is plotted as a function of the aspect ratio. For 
Xo < 2.4 the glass transition is of type B whereas for Xq > 2.5 
it is of type B' (see text). 



FIG. 5. For X,, = 2.5 and (f> = 0.593, the 

smallest aspect ratio where the type-B' transition occurs 
first the q-dependence of the non-ergodicity parameter 
Fii,(q,m) is plotted together with the static structure fac- 
tor Sw {q, m) and the normalized non-ergodicity parameter 
f[i,{q,m) = Fli,{q,m)/Siii{q,m). In (a) this is done for 
I = I' = 2 and m = whereas (b) shows the same quanti- 
ties for the center of mass I — I' = m = 0. 



FIG. 6. The formation of the 0c glass transition is 
illustrated. Within each domain of diameter ^ the system 
shows liquid crystalline order whereas for I 3> C there are 
randomly frozen orientational correlations. 
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